Importance of Granular Structure in the Initial Conditions for the Elliptic Flow 
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We show effects of granular structure of the initial conditions (IC) of hydrodynamic description of 
high-energy nucleus-nucleus collisions on some observables, especially on the elliptic-flow parameter 
112. Such a structure enhances production of isotropicaUy distributed high-pT particles, making «2 
smaller there. Also, it reduces V2 in the forward and backward regions where the global matter 
density is smaller, so where such effects become more efflcacious. 
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I. INTRODUCTION 

It is by now widely accepted that hydrodynamics is 
a successful approach for describing the collective flow 
in high-energy nuclear collisions. The basic assumption 
in hydrodynamical models is the local thermal equilib- 
rium. It is assumed that, after a complex process in- 
volving microscopic collisions of nuclear constituents, at 
a certain early instant a hot and dense matter is formed, 
which would be in local thermal equilibrium. Usually, 
this state is characterized by some initial conditions (IC), 
parametrized as smooth distributions of thermodynamic 
quantities and four-velocity. After this instant, the sys- 
tem would evolve hydrodynamically, following the well 
known set of differential equations. 

However, since our systems are not large, important 
event-by-event fluctuations are expected in real collisions. 
Concerning this question, fluctuation in the IC deserves 
a special consideration. In the past few years, we have 
studied several effects caused by such fluctuating IC on 
observables, by using a computational code especially de- 
veloped for this purpose, which we call NeXSPheRIO 
[J, i, i, i, i. In particular, in [H we have studied the 
fluctuations of the so called elliptic flow parameter V2, 
showing that they are quite large, which has been con- 
firmed by recent data [g, B- Our more recent computa- 
tions gave similar results [8]. 

An important point with regard to such IC for high- 
energy nuclear collisions is that they are not only event- 
by-event fluctuating, but also are strongly inhomogc- 
neous in space. Since the incident nuclei are not smooth 
objects, if the thermalization is verified at very early 
time as usually assumed in hydrodynamic approach, they 
could not be smooth but should have granular structure. 
In the present work, we focus our attention mainly to 
such a granular structure of fiuctuating IC and try to 
show important effects on some observables, especially 
on the elliptic-flow parameter V2. 

In what follows, we will first give a brief discussion on 
what is expected if such hot blobs are produced. This will 
be done in the next Section. We shall then describe in 
Section [TTIl our main tool, the NeXSPheRIO code. We 



show the results of computations in Section HVl first on 
transverse- momentum (pt) spectra, and then on V2 as 
function both of pseudo-rapidity 77 and of px . Finally, 
main conclusions are drawn in SeclVl 



II. 



WHAT IS EXPECTED FROM THE HOT 
BLOBS? 
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FIG. 1: Pictorial representation of the energy density distri- 
bution in a fluctuating IC. 

What do we expect if the IC present granular struc- 
ture as depicted in Fig.[T]? Because of high concentration 
of energy in point-like regions, we imagine that initially 
each blob would suffer a violent explosion and, because 
of their small size, isotropicaUy. If one of such blobs is 
deep inside the hot matter, this initial motion is quickly 
absorbed by the surrounding medium, so would not re- 
sult in any observable effect. However, if such a blob is at 
the surface of the matter, certainly the outgoing part of 
this initial acceleration would remain producing high-pT 
particles, which would be isotropicaUy distributed in the 
momentum space. 

Thus, first we expect that high-pT part of the px spec- 
tra is enhanced when fluctuating IC are used in our com- 
putations, in comparison with the results with averaged 
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(smooth) IC. In the second place, we expect that the 
anisotropic-flow coefficient V2 suffers reduction as we go 
to high-pT region, due to the additional high-pr isotropic 
component included now. As for the 77 dependence of V2 , 
we know that the average matter density decreases as 
\r]\ increases as reflected in the i] distribution of charged 
particles, so when such a blob is formed in the large-|ry| 
regions, its effects appear more clearly. Therefore, we 
expect considerable reduction of V2 in those regions. 

Although not discussed here, the granular structure we 
are considering certainly affect the so-called HBT radii. 
This question has been discussed in a previous publica- 
tion 

III. NEXSPHERIO CODE 

Our fundamental tool for the present study is called 
NeXSPheRIO. It is a junction of two computational 
codes: NeXus and SPheRIO. The NeXus code [9] is 
used to compute the IC: T^" and on some ini- 
tial hypersurface. It is a microscopic model based on 
the Regge-Gribov theory and the main advantage for 
our purpose is that, once a pair of incident nuclei or 
hadrons and their incident energy are chosen, it can 
produce, in the event-by-event basis, detailed space dis- 
tributions of energy-momentum tensor, baryon-number, 
strangeness and charge densities, at a given initial time 
r = \/W—z^ ~ 1 fm. Remark that, when we use a micro- 
scopic model to create a set of IC for hydrodynamics, the 
generated energy-momentum tensor does not necessarily 
correspond to that of local equilibrium, so we need to 
transform it to that of the equilibrated matter, adopting 
some procedure as described in detail in Ref. [lO]. 

We show in Fig. [2] an example of such a fluctuating 
event, produced by NeXus event generator, for central 
Au -I- Au collision at ISOAGeV, compared with an av- 
erage over 30 events. As can be seen, the energy-density 
distribution for a single event (left), at the mid-rapidity 
plane, presents several blobs of high-density matter, 
whereas in the averaged IC (right) the distribution is 
smoothed out, even though the number of events is only 
30. The latter would corresponds to the usually adopted 
smooth and symmetrical IC in many hydrodynamic cal- 
culations. The bumpy event structure, as exhibited in 
Fig. [21 was also shown in calculations with HIJING [ll| 
and other event generators. As already observed there 
and studied in 111, 3, i, H H, this bumpy structure 
gives important consequences in the observables. As for 
the velocity distribution, it gives essentially zero trans- 
verse velocity and longitudinal component close to the 
boost-invariant one. 

Solving the hydrodynamic equations for events, so ir- 
regular as the one shown in Fig [21 requires a special care. 
The SPheRIO code is well suited to computing the hy- 
drodynamic evolution of such systems. It is based on 
Smoothed Particle Hydrodynamics (SPH), an algorithm 
originally developed in astrophysics [1^ and adapted to 
relativistic heavy ion collisions [l3| . It parametrizes the 
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FIG. 2: Examples of initial conditions for central Au-|-Au 
collisions given by NeXus at mid-rapidity plane. Tiie energy 
density is plotted in units of GeV/fm^. Left: one random 
event. Right: average over 30 random events (corresponding 
to the smooth initial conditions in the usual hydro approach) . 

flow in terms of discrete Lagrangian coordinates attached 
to small volumes (called "particles" ) with some conserved 
quantities. Its main advantage is that any geometry in 
the initial conditions can be incorporated and this is done 
giving a desired precision. 

Now, we have to specify some equation of state (EoS) 
describing the locally equilibrated matter. Here, in ac- 
cordance with Ref. f4!| , we will adopt a phenomenological 
implementation of EoS, giving a critical end point in the 
quark-gluon plasma - hadr qn g as transition line, as sug- 
gested by the lattice QCD [14| . 

We shall neglect in this paper any dissipative effects 
and also assume the usual sudden freezeout at a con- 
stant temperature. As for the conserved quantities, be- 
sides the energy, momentum and entropy, we consider 
just the baryon number. Although very simplifled, we 
believe that the main outcomes of the present study will 
remain valid in more detailed description. 

In computing several observables, we perform in the 
present work two sets of computations: i) First, we av- 
erage over random NeXus events, obtaining smooth IC, 
which are used to compute the observables by using the 
SPheRIO code. This is similar to the usual hydro calcu- 
lation, ii) In the second set, NeXSPheRIO is run many 
times and an average over final results is performed. This 
mimics experimental conditions more closely. Remark 
that in the latter the granular structure of IC, mentioned 
above, are being explicitly included whereas in the former 
not. 

Having depicted our tool, let us now explain how we fix 
the parameters of the model and compute the observables 
of our interest. Certainly any model to be considered 
as such should reproduce the most fundamental, global 
quantities involving the class of phenomena for which it is 
proposed. So, we begin by fixing the initial conditions so 
as to reproduce properly the (pseudo-) rapidity distribu- 
tions of charged particles in each centrality window. This 
is done by applying an ry-dependent factor ~ 1 to the ini- 
tial energy density distribution of all the events of each 
centrality class, produced by NeXus. Next, we would like 
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to correctly reproduce the transverse-momentum spectra 
of charged particles, which can be achieved by choosing 
an appropriate centrality dependent freezeout tempera- 
ture, Tfo ■ 

IV. RESULTS 

A. Charged-particle spectra 

We show in Fig. [3], our results for the charged-particle 
spectra in two different centrality windows as indicated, 
computed as explained in Section IIIII It is clearly seen 
that, as one goes from the smooth averaged IC to bumpy 
fluctuating IC, the high-p^ component of the spectra 
increases as expected, making them more concave and 
closer to the data. 
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FIG. 3: Charged-particle pr distributions, in two different 
centrality windows, computed in two different ways as ex- 
plained in the text. The solid lines indicate results for event- 
by-event fluctuating IC, whereas the dotted lines the ones for 
the averaged IC. The data points are also plotted for compar- 
ison [l5| . 



B. pt dependence of (172) 

In Fig.|4], we show our results for the pt dependence 
of {V2) in the centrality window and 77 interval as indi- 
cated. Here, since our purpose is to show clearly the 
effects of granular IC, we plot the results obtained with 



the same PHOBOS hit-based method for the both curves, 
but without any correction for event-plane fluctuations. 
While the result with smooth averaged IC shows contin- 
uous increase of (^2) with px , deviating largely from the 
data points, as happens in the orthodox hydro compu- 
tations, the introduction of spiky IC makes V2 smaller 
at high Pt as expected, approaching the curve to the 
data points. The correction mentioned above shifts up- 
ward the curve, corresponding to the fluctuating IC, but 
without modifying its shape. The fact that the averaged 
smooth IC lead to rising V2{pt) for large pt and not flat- 
tening is usually interpreted as a breakdown of hydrody- 
namic model. We suggest it could be in part related to 
the granular structure of the IC. 
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FIG. 4: Pt dependence of (112) in the centrality window and 
77 interval as indicated, compared with data [l6l ]. The solid 
line indicates result for fluctuating IC, whereas the dotted one 
that for the averaged IC. The curves are averages over PHO- 
BOS centrality sub-intervals with freeze-out temperatures as 
indicated. 



C. r\ dependence of (112) 

In Fig. El, we show our results for the ry dependence of 
(f 2) in three different centrality windows as indicated. As 
in Ref. we calculated {v-i) with respect to the event 
plane as done experimentally. The event plane has been 
determined here by all the charged particles (the results 
with the PHOBOS hit-based method with correction are 
almost identical). For averaged smooth IC, in agree- 
ment with the usual hydro computations [13, [isl. Il9||. 
{v-i) exhibits shoulders in high-jTyj regions. When ffuctu- 
ating spiky IC are used, these shoulders are considerably 
weakened as expected. The results, which combine this 
and the complementary effect of increase of {v-i) with the 
event-plane fluctuations, are now very close to the data. 

Another important ingredient of a hydrodynamic 
model is the decoupling procedure. We are studying the 
effect of continuous emission [20| . Probably it makes the 
curve of V2 at high-p^ even flatter and the jy-distribution 
narrower [4|. This is clear, also in Refs. [H, [H, [2]| . 
where they use the usual smooth IC and, by describing 
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FIG. 5: ri dependence of (U2) for three centrality windows. 
The soUd Unes indicate results for event-by-event fluctuating 
IC, whereas the dotted hnes the ones for the averaged IC. The 
data points are also plotted for comparison [l^ . T/o has been 
taken as indicated in FIG|31 

the decoupling with a microscopic transport model, could 
improve V2 in this direction. The cascade model in the 
hadronic stage produces similar effects as the continuous 
emission. 



IC, the expectation being that the granularity will not 
manifest itself. In this paper, we argue that this may not 
be true qualitatively and also quantitatively. 

The main conclusions of the present study are: 

1. Granular structure in the IC produces more high- 
Pt particles, distributed isotropically. 

2. Granular structure in the IC reduces the elliptic 
flow, because of the isotropic flow it generates. 

3. As function of px, the mechanism becomes more 
effective as px increases, because those high-density 
blobs cause violent expansion, producing high-p^ 
particles. 

4. This effect is enhanced where the average matter 
density is small. So, it decreases V2 in the large 
pseudorapidity regions. 

5. NeXSPheRIO, with fluctuating and spiky IC, re- 
produce approximately the data both of the pr- 
distribution and jy-distribution of V2 for different 
centrality windows. 



V. CONCLUSIONS 

While it is believed that a granular structure exists in 
the IC, hydro simulations are usually done with smooth 
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